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The high energy neutrino cross section is a crucial ingredient in the calculation of the event rate 
in high energy neutrino telescopes. Currently there are several approaches which predict different 
behaviours for its magnitude for ultrahigh energies. In this paper we present a comparison between 
the predictions based on linear DGLAP dynamics, non-linear QCD and in the imposition of a 
Froissart-like behaviour at high energies. In particular, we update the predictions based on the 
Color Glass Condensate, presenting for the first time the results for ovjv using the solution of the 
running coupling Balitsky-Kovchegov equation. Our results demonstrate that the current theoretical 
uncertainty for the neutrino-nucleon cross section reaches a factor three for neutrinos energies around 



PACS numbers: 12.38.-t; 13.10.+q; 13.15.+g 



I. INTRODUCTION 

The study of ultra high energy (UHE) cosmic neutrinos 
provides an opportunity for study particle physics beyond 
the reach of the LHC In particular, the Pierre Auger 
Observatory (PAO) is sensitive to neutrinos of energy 
> 10 8 GeV [2]. A crucial ingredient in the calculation 
of attenuation of neutrinos traversing the Earth and the 
event rate in high energy neutrino telescopes is the high 
energy neutrino-nucleon cross section, which provides a 
probe of Quantum Chromodynamics (QCD) in the kine- 
matic region of very small values of Bjorken-x. The typ- 
ical x value probed is x re m,y V -/2m^rE I/ , which implies 
that for E v re 10 s - 10 10 GeV one have x re 1(T 4 - 1(T 6 
at Q 2 re 10 4 GeV 2 . This kinematical range was not ex- 
plored by the HERA measurements of the structure func- 
tions n. 

The description of QCD dynamics in the high energy 
limit still is a subject of intense debate (For a recent 
review see e.g. Ref. Q). Theoretically, at high en- 
ergies (small Bjorken-x) one expects the transition of 
the regime described by the linear dynamics, where only 
the parton emissions are considered, to a new regime 
where the physical process of recombination of partons 
becomes important in the parton cascade and the evo- 
lution is given by a non-linear evolution equation. This 
regime is characterized by the limitation on the maxi- 
mum phase-space parton density that can be reached in 
the hadron wavefunction (parton saturation), with the 
transition being specified by a typical scale, which is 
energy dependent and is called saturation scale Q s 
Moreover, the growth of the parton distribution is ex- 
pected to saturate, forming a Color Glass Condensate 
(CGC), whose evolution with energy is described by an 
infinite hierarchyof coupled equations for the correlators 
of Wilson lines [!, Hj] . In the mean field approximation, 
the first equation of this hierarchy decouples and boils 
down to a single non-linear integro-differential equation: 
the Balitsky-Kovchegov (BK) equation [U, 0- Experi- 
mentally, possible signals of parton saturation have al- 
ready been observed both in ep deep inelastic scattering 



at HERA and in deuteron-gold collisions at RHIC (See, 
e.g. Ref. iil). 

Existing estimates of the neutrino nucleon cross sec- 
tions with structure functions constrained by HERA data 
are based on linear dynamics (DGLAP or an unified 
DGLAP/BFKL evolution) or phenomenological models 
that mimics the expected behaviour predicted by the 
non-linear QCD dynamics. In particular, the neutrino - 
nucleon cross section was originally calculated at leading 
order in Ref. [To| , with the resulting parametrization be- 
ing a benchmark for the evaluation of sensitivies of UHE 
cosmic neutrinos. In Refs. [ll|, Ell a next-to-leading or- 
der analysis was performed, and the uncertainties on high 
energy <7„at which are compatible with the conventional 
DGLAP formalism [13[ was estimated. Moreover, in Ref. 
14] it was estimated considering an analytical solution 
of the DGLAP equation, valid at twist- 2 and small-x. 
These approaches imply a power-like increase of the cross 
section at ultrahigh energies which is directly associated 
to the DGLAP solution at small-x. In contrast, in Refs. 
[l5l HH ] the HERA data were successfully fitted assuming 
that the proton structure function saturates the Froissart 
bound, which imp lies F"^ oc ln 2 (l/x). On the other hand, 
in Refs. |17H20l | the contribution of non-linear effects 
was estimated considering the leading order solution of 
the BK equation and distinct phenomenological models 
based on saturation physics [H, [l9j]. These non-linear 
approaches were improved recently with the solution of 
the BK equation including running coupling corrections 
and the successful description of the HERA data. This 
fact motivates a revision of previous estimates. In partic- 
ular, we present for the first time the predictions for the 
neutrino-nucleon cross section obtained using as input 
the numerical solution of the Balitsky-Kovchegov equa- 
tion 0, Q including running coupling corrections |2l| - 
[23|. Moreover, we present the predictions of some recent 
phenomenological parametrizations based on saturation 
physics, which provide an economical description of a 
wide range of data with a few parameters. Our main mo- 
tivation is to provide an update on the neutrino-nucleon 
cross section in the literature and estimate the theoretical 
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uncertainty in its behaviour at ultra high energies. 

This paper is organized as follows. In next section (Sec- 
tion [TTJ we present a brief review of the neutrino-nucleon 
deep inelastic scattering (DIS) cross section, introducing 
the main formulae. In Section Mil we discuss the QCD 
dynamics and the models used in the calculations. In 
Section IIVI we compare the predictions of the different 
models. Finally, in Section |V] we summarize our main 
results and conclusions. 



II. THE NEUTRINO - NUCLEON DIS CROSS 
SECTION AT HIGH ENERGIES 

Deep inelastic neutrino scattering is described in terms 
of charged current (CC) and neutral current (NC) inter- 
actions, which proceed through and Z° exchanges, 
respectively. The total cross sections are given by [24| 



CC, NC r jTi \ 
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where E v is the neutrino energy, s — 2ME l/ with M the 
nucleon mass, y = Q 2 /(xs) and Ql, lin is the minimum 
value of Q 2 which is introduced in order to stay in the 
deep inelastic region. In what follows we assume Q 2 m i n = 
1 GeV 2 . Moreover, the differential cross section is given 
by 
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where Gp is the Fermi constant and Mj denotes the mass 
of the charged of neutral gauge boson. The calculation 
of o v n involves integrations over x and Q 2 . On the one 
hand, as the parton distributions behaves as x~ x (A > 0) 
at low x, the x integral becomes dominated by the in- 
teraction with partons of lower x. On the other hand, 
the Q 2 integral remains dominated by Q 2 values to the 
order of the electroweak boson mass squared. For Q 2 
above the integrand behaves as 1/Q 4 and quickly 
becomes irrelevant. As the neutral current (NC) interac- 
tions are subdominant, we will consider in what follows, 
for simplicity, only charged current (CC) interactions. 

In the QCD improved parton model the structure func- 
tions Fi(x, Q 2 ) are expressed in terms of the parton dis- 
tributions on the nucleon, which satisfy the DGLAP 
fl3j and/or BFKL [25[ linear dynamics. On the other 
hand, an efficient way of introducing non-linear effects 
is the description of the structure functions considering 
the color dipole approach in which the DIS to low x 
can be viewed as a result of the interaction of a color 
qq dipole which the gauge boson fluctuates [2^|. In this 
approach the F 2 ' structure function is expressed in 



terms of the transverse and longitudinal structure func- 
tions, F? C ' NC = F° C - NC + F? C > NC which are given 
by 



F, 



CC,NC 
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where r denotes the transverse size of the dipole, z is 
the longitudinal momentum fraction carried by a quark 
and ^ T '^ are proportional to the wave functions of the 
virtual charged of neutral gauge bosons corresponding to 
their transverse or longitudinal polarizations. Explicit 
expressions for ^^'^ are given, e.g., in Ref. [17| . Fur- 
thermore, a dp describes the interaction of the color dipole 
with the target. In next section we will discuss some 
models for <j dp , based on the non-linear QCD dynamics, 
which describe the current HERA data. 

A comment is in order. The HERA measurements of 
the structure function at low - x (x w 10~ 6 ) are for 
very low values of Q 2 (Q 2 -C 1 GeV 2 ), which implies 
that small x extrapolations of the parton distributions 
are necessary to estimate <j v n above E v w 10 7 GeV. As 
we will discuss in more detail in the next sections, these 
extrapolations contain significant uncertainties. 



III. HIGH ENERGY QCD DYNAMICS 

High energy neutrino-nucleon cross section accesses 
very large values of Q 2 and very small values of Bjorken x. 
In the last years, it has been calculated considering differ- 
ent theoretical approaches with structure functions con- 
strained by HERA data. It was originally calculated us- 
ing the LO DGLAP equation in Ref. OH (GQRS model). 
In Ref. [IH (C-SS model) a next-to-leading order anal- 
ysis was performed considering a global fit of the ZEUS 
data and the uncertainties on high energy <7„jv which are 
compatible with the conventional DGLAP formalism was 
estimated. Another approach based on DGLAP dynam- 
ics was proposed in [lj] (FJKPPP model), where ovjv 
was estimated using an analytic result for the DGLAP 
evolution of the structure functions, valid at twist-2 in 
the region of small-x and for a soft non-perturbative in- 
put. As the DGLAP equation can break down at low x 
because of the ln(l/x) terms which appear in the per- 
turbative series, an alternative approach which incorpo- 
rates both the ln(l/x) resummation and the complete 
LO DGLAP evolution was proposed in (27j, which we 
denote unified DGLAP-BFKL model hereafter (See also 
Ref. [H). 

These predictions are based on the linear DGLAP 
and/or BFKL equations which implies a power increase 
with the energy of the neutrino-nucleon cross section that 
eventually violate the Froissart bound. In Ref. [HI the 
UHE a V M was estimated using a phenomenological ap- 
proach (BBMT model) which considers the imposition 
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of the Froissart unitarity and analyticity constraints on 
inclusive deep-inelastic cross sections. It implies that 
F% oc ln 2 (l/a;) at very small x. In fla ] very good fits 
to data were obtained for x < 0.1 and a wide range of 
Q 2 . More recently, in [l6j], this model was updated con- 
sidering the most recent analysis of the complete ZEUS 
and HI datasets from HERA. In what follows we denote 
by BHM the corresponding predictions. 

Another alternative to estimate the ultrahigh energy 
behavior of the neutrino-nucleon cross section is to ex- 
press the structure functions in the dipole approach [Eq. 
(J3j> ] and to consider the state-of-art of the non-linear 
QCD dynamics: the Color Glass Condensate formalism. 
In this formalism, the dipole - target cross section a dp 
can be computed in the eikonal approximation, resulting 



a dp (x,r) = 2 J d 2 bM{x,r,b) 



(4) 



where Af is the dipole-target forward scattering am- 
plitude for a given impact parameter b which encodes 
all the information about the hadronic scattering, and 
thus about the non-linear and quantum effects in the 
hadron wave function. It is useful to assume that the 
impact parameter dependence of J\f can be factorized as 
Af(x,r,b) =N(x,r)S(b), so that a dp (x,r) = a Af(x,r), 
with (To being a free parameter related to the non- 
perturbative QCD physics. The Balitsky-JIMWLK hier- 
archy describes the energy evolution of the dipole-target 
scattering amplitude Af(x, r). In the mean field approx- 
imation, the first equation of this hierarchy decouples 
and boils down to the Balitsky-Kovchegov (BK) equa- 
tion 0,0. 

In the last years the next-to-leading order corrections 
to the BK equation were calculated [21H231 ] through the 
ressumation of a s Nf contributions to all orders, where 
Nf is the number of flavors. Such calculation allows one 
to estimate the soft gluon emission and running coupling 
corrections to the evolution kernel. The authors have 
verified that the dominant contributions come from the 
running coupling corrections, which allow us to deter- 
mine the scale of the running coupling in the kernel. The 
solution of the im pro ved BK equation was studied in de- 
tail in Refs. [22I. |29||. Basically, one has that the run- 
ning of the coupling reduces the speed of the evolution 
to values compatible with experimental data, with the ge- 
ometric scaling regime being reached only at ultra-high 
energies. In [30] a global analysis of the small x data 
for the proton structure function using the improved BK 
equation was performed (See also Ref. [3l[). In contrast 
to the BK equation at leading logarithmic a s ln(l/x) ap- 
proximation, which fails to describe the HERA data, the 
inclusion of running coupling effects in the evolution ren- 
ders the BK equation compatible with them (See also 

[HH). 

The dipole-target cross section can also be calcu- 
lated considering phenomenological parametrizations for 
Af(x, r) based on saturation physics, which provide an 
economical description of a wide range of data with a 
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FIG. 1: (Color online) The neutrino nucleon CC cross section 
as a function of the neutrino energy, E v . Comparison between 
the linear and Froissart-inspired approaches. 



few parameters. Several models for the forward dipole 
cross section have been used in the literature in order to 
fit the HERA and RHIC data @,|HH4f|. In general, the 
dipole scattering amplitude is modelled in the coordinate 
space in terms of a simple Glaubcr-like formula as follows 



N(x, r) = 1 — exp 



(5) 



where 7 is the anomalous dimension of the target gluon 
distribution. The main difference among the distinct 
phenomenological models comes from the predicted be- 
haviour for the anomalous dimension, which determines 
the transition from the non-linear to the extended ge- 
ometric scaling regimes, as well as from the extended 
geometric scaling to the DGLAP regime (See e.g. jH). 
The current models in the literature consider the general 
form 7 = 7 S + A7, where j s is the anomalous dimen- 
sion at the saturation scale and A7 mimics the onset of 
the geometric scaling region and DGLAP regime. One of 
the basic differences between these models is associated 
to the behaviour predicted for A7. While the models 
proposed in Refs. [i^ - li^ ] assume that A7 depends on 
terms which violate the geometric scaling, i.e. depends 
separately on r and rapidity Y = ln(l/x), the model 
proposed in Ref. [45[ (BUW model) consider that it is 
a function of rQ s . In particular, these authors demon- 
strated that the RHIC data for hadron production in dAu 
collisions for all rapidities are compatible with geometric 
scaling and that geometric scaling violations are not ob- 
served at RHIC [45j]. In contrast, the IIM analysis [42[ 
implies that a substantial amount of geometric scaling 
violations is needed in order to accurately describe the 
ep HERA experimental data. In our analysis we will esti- 
mate a V N considering the rcBK solution and the IIM and 
BUW models. Moreover, we will consider the improved 
version of the IIM model proposed in Ref. [46| (denoted 
hereafter IIM-S model), which includes the heavy quark 
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FIG. 2: (Color online) The neutrino nucleon CC cross section 
as a function of the neutrino energy, E v . Comparison between 
the non-linear approaches. 
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FIG. 3: (Color online) The neutrino nucleon CC cross section 
as a function of the neutrino energy, E v . Comparison between 
the linear and non-linear approaches. 



effects on the saturation in the fit of the HERA data. The 
main differences of the IIM-S model in comparison to the 
IIM one are the larger value of the anomalous dimension 
and the smaller value of the exponent which determines 
the energy growth of the saturation scale. Moreover, the 
IIM-S model considers a newer HI and ZEUS datasets in 
the fit, in contrast to the IIM one which only considers 
the ZEUS data. 



IV. RESULTS 

In this section we present a comparison between the 
predictions of the linear approaches (GQRS, C-SS and 
FJKPPP), the Froissart-inspired models (BBMT and 
BHM) and non-linear approaches (rcBK, IIM, IIM-S and 
BUW). In Fig. [T]the energy dependence of the neutrino 
nucleon CC cross section predicted by the linear and 
Froissart-inspired approaches are compared. As expected 
from the solution of the DGLAP equation at small- a;, 
the GQRS, C-SS and FJKPPP models predict a strong 
increase of the cross section at ultrahigh energies. Al- 
though these approaches agree at low energies, where the 
behavior of the parton distributions are constrained by 
the HERA data, they differ by a factor 1.25 at E v = 10 12 
GeV. We have that the C-SS prediction, which comes 
from a global fit of the ZEUS data, can be considered as 
a lower bound for the linear predictions. On the other 
hand, the FJKPPP prediction, which considers an ana- 
lytical solution of the DGLAP equation at small-x, im- 
plies a stronger increase with the energy similar to the 
GQRS one, largely used to estimate the event rates in 
neutrinos telescopes. It is important to emphasize that 
GQRS parameterization was obtained using a restrict set 
of experimental data and DGLAP evolution at leading 
order. In contrast, the BBMT and BHM approaches, 
which assume a Froissart-like behavior for the structure 



functions at small- a; [F2 oc ln 2 (l/a;)], predict at ultra high 
energies a cross section smaller than the FJKPPP one by 
a factor « 3. We have that the BBMT and BHM predic- 
tions, which comes from a fit of the structure functions 
using different datasets, are very similar. 

In Fig. [5] we present the predictions of the non-linear 
approaches for the energy dependence of the neutrino nu- 
cleon CC cross section. For comparison, we also includes 
the GBW model [35[ which was used in previous stud- 
ies of the saturation effects [l?], [01 ■ We have that the 
IIM-S prediction is very similar to the GBW one, being 
a lower bound for the non-linear predictions at ultrahigh 
energies. In contrast, the rcBK prediction can be consid- 
ered an upper bound. These predictions differ by a factor 
ss 3 at E v = 10 12 GeV. The large difference between the 
IIM and IIM-S predictions is directly associated to the 
distinct energy dependence for the saturation scale and 
the treatment of heavy quarks proposed by these models. 
It implies a different normalization for the dipole-target 
cross section and, consequently, for the neutrino-nucleon 
CC cross section. The IIM, BUW and rcBK predictions 
are similar at low energies, but differ by a factor 1.8 at 
ultra high energies. It is important to emphasize that 
BUW and rcBK models successfully describe the current 
RHIC and HERA data 0, S3, 

Finally, in Fig. [3] we present a comparison between 
the predictions of linear and non-linear approaches. For 
comparison we also include the prediction obtained by 
the unified DGLAP-BFKL approach < 27J, which is sim- 
ilar to the FJKPPP one. We have that the FJKPPP 
and rcBK predictions are similar for E v = 10 11 GeV and 
differ by ss 15% at E v = 10 12 GeV. In contrast, the 
rcBK and BBMT differ by a factor ss 3 at E v = 10 11 
GeV. The theoretical uncertainty increases for a factor 
ss 5.5 when we compare the FJKPPP and BBMT pre- 
dictions for E v — 10 13 GeV. Our results demonstrate 
that the determination of <j v n can be useful to contrain 
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the underlying QCD dynamics. In principle, this cross 
section could be constrained at high energies by study- 
ing the ratio between quasi-horizontal deeply penetrating 
air showers and Earth-skimming tau showers [ill ]. 

V. SUMMARY 

Detection of UHE neutrinos may shed light on the ob- 
servation of air showers events with energies in excess of 
10 11 GeV, reveal aspects of new physics as well as of the 
QCD dynamics at high energies. One of the main ingredi- 
ents for estimating event rates in neutrino telescopes (e.g. 
ICECUBE) and cosmic ray observatories (e.g. AUGER) 
is the neutrino - nucleon cross section. In this paper we 
examined to what extent the cross section is sensitive to 
the presence of new dynamical effects in the QCD evolu- 
tion. We compare the predictions of several approaches 
based on different assumptions for the QCD dynamics. 
In particular, we have compared the more recent predic- 



tions based on the NLO DGLAP evolution equation with 
those from the CGC physics obtained using the running 
coupling BK solution or phenomcnological models. Our 
results demonstrate that the current theoretical uncer- 
tainty for the neutrino-nucleon cross section reaches a 
factor three for neutrino energies around 10 11 GeV and 
increases to 5.5 for E v — 10 13 GeV. 

A final comment is in order. In this paper we esti- 
mated the range of possible values for er„jv at large en- 
ergies within the Standard Model. It makes possible to 
search for enhancements in the neutrino-nucleon cross 
section due to physics beyond the perturbative SM (See, 
e.g. 0). 
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